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The dispersion-managed (DM) optical system with step- wise periodical variation of dispersion is 
studied in the framework of path-averaged Gabitov-Turitsyn equation. The soliton solution is ob- 
tained by iterating the path-averaged equation. The dependence of soliton parameters on dispersion 
map strength is investigated together with the oscillating tails of soliton. 
OCIS codes: 060.2330, 060.5530, 060.4370, 190.5530, 260.2030. 



Dispersion-managed (DM) systemB, which is the sys- 
tem with periodical variation of dispersion along opti- 
cal fiber, is one of the key point of current developraenjb 
of ultrafast high-bit-rate optical communication linesErQ. 
The main factor limiting the bit-rate capacity is pulse- 
broadening due to the chromatic dispersion of the op- 
tical fiber which can be overcome by periodical chang- 
ing of sigh of fiber's dispersion to create very low (or 
even zero) path-average dispersion. Lossless propaga- 
tion of optical pulse in DM fiber is described by nonlin- 
ear Schrodinger equation (NLS) with periodically vary- 
ing dispersion d{z): 



d{z)utt 



-u = 0, 



(1) 



where u is the envelope of optical pulse, z is the prop- 
agation distance and all quantities are made dimen- 
sionless. Consider two-step periodic dispersion map: 
d{z) = (d) + d{z), where d{z) = di for < z + < L/2 
and d{z) — —di for L/2 < z + nL < L, (d) is the path- 
averaged dispersion, di is the amplitude of dispersion 
variation, L is a dispersion period and n is an arbitrary 
integer number. Eq. (^ also describes pulse propagation 
in fiber with losses compensated by periodically placed 
amplifiers if the distance between amplifiers is much less 
than L. 

Provided a characteristic nonlinear length Zni of the 
pulse is large: Zni ^ L, where Zni = and p is a 

typical pulse amplitude, the eq^, (|l]) can be reduced to 
path-averaged Gabitov-TuritsynQ model: 
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good supported by muperical simulationsE 
to t-space one can geta: 



Returning 



ii'z + {d)iptt 
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Cii — )^{h + t)^{t2 + t) X 

s 

tp* {ti + t2 + t)dtidt2 = 0, (3) 



where Ci{x) — ^^^dx Lnote difference in definition of 
Ci{x) in comparison withal). It was found numericallycl 
that Gaussian ansatz 



A-Gauss = P exp I 
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(4) 



ip{uJ2 + {<-^Jl + UJ2 + uj)dujiduj2 = 0, (2) 

where s = diL/2 is a dispersion map strength, ip = 

„ iuj'^ i d[z')dz' ^ , . /.OC ; / \ 7 i# 7 • -n 

ue •'^1^ and ip\Lo) = J _^ip[t)e dt is Courier 

component of ip. Gabitov-Turitsyn model is based on 
the assumption of slow variation of ^ as a function of z 
on scales of the order of DM period L. This model is 



where p, P are real constants, is a rather good approxi- 
mation for the soliton solution u — A{t)e^^^ {A is real) 
of eq. (|^) for small and moderate values of t at space 
points z = L/4 + nL. Respectively (Q) is also a good 
approximation for the soliton solution ip — A{t)e^'^^ of 
(^. On soliton solution eq. (^) takes the form: 

-XA+ {d)Au = ^ / CiC-^)Aih +t)x 

A{t2 + t)A{ti + t2 + t)dhdt2, (5) 

Based on this observation the soliton solution of (|l|) 
was approximated inll30 by Gaussian exponent multi- 
plied by a sum consisting of a finite number of Hermite 
polynomials. Because infinite set of these polynomials 
is complete in the space of square-integrable function, ^ 
can be expanded over this set. But for large enough 
t such approximatioiU|Can not be effective because of 
highly oscillating tailsBllj of A for t ^ oo which requires 
to take into account very large number terms of expan- 
sion to get reasonable approximation. The oscillating 
tails of A are of big importance because they are respon- 
sible for interaction of a sequence of solitona launched in 
optical fiber which limits bit-rate capacityO. 

Here the convergent set of approximate solutions of (^) 
is found by iterating that eq. for (d) — 0. Zero iteration 
is given by (0). — AA'"^ is obtained by substitution 
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of into right hand side of (|) for n = 1, 2, . . . The 

oscillating tails appear at first iteration already. Simul- 
taneously one can get the dependence of [3 in eq. (|^) on 
dispersion map strength s. 

Rewriting the kernel of (p|) via parametric integral: 



sinsa;iCi;2 1 
SLO1LO2 2s 



(6) 



and using as zero iteration allow to make integra- 
tion over time variables ti, t2 in (p|) explicitly and get 
expressions for first A^-'^' , second A'^^ , etc. iterations: 



exp I 



VI - 2zs' + 3s'2 
exp(-gt2) ^ 
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ds dsids2ds3, (7) 



where s 



f3s, integration limits for all variables are 
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(&1 + 62 + &3)s'2), 



1,2,3, 5 = 4(46i62&3 - 416162s' 

q = (4(6162 + 6163 + 6263) + 4is'63 + s'2)/g, and p^") is 
determined by the condition A'^"^|i=o — V for all n. Note 
that instead of keeping constant for all itera- 

tions it is also possible to fix the integral Pq = / |^^"^ 
if one aims to get soliton solution with some definite 
value of Pq. In principle any values of A, /3 can be chosen 
for zero iteration but the most useful is to choose them 
in such a way to get A^^^ as close as possible to A^^^ in 
order to get faster convergence of iterations to true solu- 
tion of eq. (^. It is convenient to make series expansion 
Ad) in powers of which allows to integrate all terms 
of this expansion explicitly. In particular 



p'-^-' = \fiy I {arcsinh—— V c.c. 
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+ c.c. 



where c.c. means complex conjugation. Equating terms 
proportional in series expansion of A^^) and (^) one 
gets from (0) transcendental equation for s which gives 



S = (3s^ 2.393 . . . 
and using (|^) with the same accuracy 
A = p2 . 0.482 . . . 



(9) 



(10) 



2.5, (d) = 0, t = i/Vs and P, A are given by (|), (|T^). 
It is enough to consider only positive values of t be- 
cause A{t) is an even function of t. Note that in case 
(d) = eq. (|^) is invariant under two scaling trans- 



mit, s 



A{t) 



A{t/Si) and 



forms: t 

A (5f A, A S2A, which allows to extend results of 
Fig.l to all values of system parameters p, s (of course it 
has a sense to consider only such values of p, s for which 
condition Zni S> L holds). 

It is seen in Fig. 1 that the convergence of iteration 
procedure to soliton solution is very fast for small t and 
becomes slower for larger t. For \i\ ^ 1.3 the renormal- 
ization of A^^) by next iteration is very small. In par- 
ticular for aU n = 2, 3, . . . (Aj"V^^") L^o differs from /3 
determined by (^ for less that 2%. For \i\ ~ 2.3 the good 
approximation is given by second iteration, for \i\ ^ 3.2 
by third etc. The sequence of iteration can be inter- 
preted as evolution of A along some artificial coordinate 
z and it can indicate that the convergence of the solution 
of (|^) during evolution along coordinate z to soliton so- 
lution (^) should also be slower for larger t what is really 
observed in nutaerical experiments if tp\z=o is chosen in 
Gaussian formllj. 

Thus, provided the form of soliton is close to Gaus- 
sian for small and moderate t, the value of /3s is unt. 
versal and is given by Zakharov and Manakovfij 

proposed the other theory of DM soliton in strong DM 
limit based on the assumption that typical width of A 
distribution is much less than ^/s. In present notation it 
means /3s 3> 1. Apriori one can not exclude that soliton 
solution with /?s 3> 1 is also possible in addition to ob- 
tained here. But such solution can not be Gaussian-like 
and, according to best of my knowledge, in all numeri- 
cal experimeniiSU-far the value of /3s was of the order of 
unity (see e.g.H&E£l). Note that choice of /3s ^ 1 for zero 
iteration results in convergence of iteration sequence (0) 
to Gaussian-like solution with /3 is given by (||). It can 
indicate that solution of (^ for /3s 3> 1 also converges for 
large z to that Gaussian-like solution. One can not also 
exclude the existence of soliton with /3s ;» 1 for (d) ^ 
which is outside the scope of the present paper. 

Let us try to explore in more details in what sense soli- 
ton solution A is close to Gaussian. From (|^) , (||) and (^) 
one can get series expansion (the approximate numeri- 
cal value of s is used here to avoid writing a cumbersome 
expression for every term of expansion) : 

^(1) exp (-^) /p ~ 1 - 8.211 • 10-2 ? - 2.654 • 10"^ P 



4.383 • 10"^?* - 3.946 • lO"'' t 



4 no 



(11) 



Fig.l shows the results of first four iterations obtained 
numerically for p — 1.7, di = 500, L = 0.01, s = 



One can conclude from this expansion that A can be 
written as multiplication of Gaussian exponent (0) on 
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a slow function which changes significantly on scales 
P ~ 3.5. For this and larger scales there is no simi- 
larity between soliton and (^). Thus DM soliton is close 
to Gaussian exponent only with accuracy up to numer- 
ically " small parameter" ~ f /4 (because scale of Gaus- 
sian exponent is (31^/2 = ^ 1) and there is no 
really small parameter describing this closeness. Never- 
theless appearance of numerical "small parameter" 1/4 
explains success of eptBiaBsion of DM soliton solution in 
Hermite polynomialsll^Hiil which is basically equivalent to 
expansion (|ll| ) in Taylor series. 

For larger scales \t\ ^ 1.3 the correct presentation of 
DM soliton requires a lot of terms of series expansion. 
Instead one can get asymptotical behavior from (0). E.g. 
the asymptotic of the first iteration is given by: 



n(l) 



4i(s - i) + 



rexp ( 
(31 
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t2 



+ c.c. + 0{^) (12) 



Fig. 2 shows that eq. ( p^ ) (curve 2) is a relatively good 
approximation of A^^"> (curve 1) for \t\ ~ 1.6. But of 
course to improve accuracy of soliton solution approxi- 
mation it is necessary to get asymptotic of next itera- 
tions. Another way is to use the series expansion ( pT| ) 
as zero iteration instead of (|^). E.g. keeping first four 
terms of this expansion for zero iteration one can get 
parametric integral similar to (^ which gives good ap- 
proximation of A^^^f in (^ for \i\ ~ 3.5. Detail consider- 
ation of this approximation is outside the scope of the 
present paper. 
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Figure captions: 



Fig.l. Time-dependence of first, second, third and fourth 
iteration (curves 1,2,3,4 respectively). 
Fig. 2. First iteration A^^'>{t) (curve 1) versus eq. ( p^ ) 
(curve 2). . 
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